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ABSTRACT 

Using cosmological hydrodynamic simulations we measure the mean transmitted flux 
in the Lya forest for quasar sightlines that pass near a foreground quasar. We find 
that the trend of absorption with pixel-quasar separation distance can be fitted using 
a simple power law form including the usual correlation function parameters tq and 7 
so that ((F(r)) = ex P(~ T off(l + (^)~ 7 )))- From the simulations we find the rela- 
tion between tq and quasar mass and formulate this as a way to estimate quasar host 
dark matter halo masses, quantifying uncertainties due to cosmological and IGM pa- 
rameters, and redshift errors. With this method, we examine data for ~ 3000 quasars 
from the Sloan Digital Sky Survey (SDSS) Data Release 3, assuming that the effect 
of ionizing radiation from quasars (the so-called transverse proximity effect) is unim- 
portant (no evidence for it is seen in the data.) We find that the best fit host halo 
mass for SDSS quasars with mean redshift z — 3 and absolute G band magnitude 
-27.5 is log lo M/M = 12.48t u ;gg. We also use the Lyman-Break Galaxy (LBG) and 
Lya forest data of Adelberger et al in a similar fashion to constrain the halo mass of 
LBGs to be log 10 M/M0=11.131q g§, a factor of ~ 20 lower than the bright quasars. 
In addition, we study the redshift distortions of the Lya forest around quasars, using 
the simulations. We use the quadrupole to monopole ratio of the quasar-Lya forest 
correlation function as a measure of the squashing effect. We find that this does not 
have a measurable dependence on halo mass, but may be useful for constraining cosmic 
geometry. 
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1 INTRODUCTION 

The Lyman-alpha forest of absorption lines seen in the spec- 
tra of quasars see (e.g., Rauch 1998 for a review) can be 
related in theories of structure formation to fluctuations in 
the matter distribution. Because the fluctuations are only 
weakly non-linear, the photoionized gas where they arise 
traces the dark matter distribution that dominates the grav- 
itational potential very well, and the Lya forest can be used 
to give us information on structure in the dark matter (e.g., 
Cen et al. 1994; Zhang, Anninos, & Norman 1995; Petit- 
jean, Miicket, & Kates 1995; Hernquist et al. 1996; Katz 
et al. 1996; Wadsley & Bond 1997; Theuns et al. 1998; Dave 
et al. 1999). We can therefore use the Lya forest absorp- 
tion around galaxies and quasars to tell us about the profile 
of dark matter around these objects. The clustering of dark 
matter around dark matter halos is expected to be related to 
their mass. This has been used by e.g., Seljak et al. (2005) 
and Mandelbaum et al. (2006), to constrain the mass of 
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galaxy halos using the matter distribution around them in- 
ferred from weak lensing observations. In the present paper, 
we propose to infer related constraints on the mass of quasar 
hosts and high redshift galaxies from the Lya forest-derived 
mass profiles around them. 

Some current mass constraints for quasar hosts and high 
redshift galaxies come from considering the amplitudes of 
their two point-correlation functions. For example, Croom 
et al. 2005 analyzed the autocorrelation function of quasars 
in 2dF data, and Myers et al. (2006) the clustering of 300,000 
photometrically classified SDSS QSO to find constraints on 
quasar host masses. 

It has been known since the work of Bajtlik et al. (1988) 
that the Lya forest in the spectrum of quasars shows evi- 
dence of decreasing absorption as the quasar's emission red- 
shift is approached. This is known as the proximity effect, 
and is believed to be due to photoionization from the quasar 
itself adding to the mean background photoionizing radia- 
tion field and decreasing the amount of neutral hydrogen in 
the vicinity of the quasar. Knowing quasar observed lumi- 
nosities, the observed decrease in Lya absorption has been 
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used to estimate the intensity of the mean cosmic back- 
ground photoionizing radiation field (see e.g., Scott et al. 
2000). A related process is the possible effect of other (fore- 
ground) quasars on the Lya forest of other quasars whose 
sightlines pass close by. This effect, known as the transverse, 
or foreground proximity effect has not so far been observed, 
and in fact there have been several non-detections of the 
effect (e.g. Schirber & Miralda-Escude 2004, Croft 2004 ). 

This could be explained if for example quasar radiation 
is beamed, so that the usual proximity effect would occur 
but off-axis sightlines would not be affected. Alternatively, 
the quasar lifetime could be short enough that light has 
not had time to travel transversely to adjacent sightlines 
(for example if they are 10 h~ Mpc away, a quasar lifetime 
10 Myr would leave them unaffected. The two studies re- 
ferred to above point instead to increased Lya absorption in 
the parts of sightlines that pass close to foreground quasars. 
This is expected, as quasars should be hosted by massive 
dark matter halos which are in overdense regions. We will 
use this fact along with theoretical predictions in the con- 
text of the Cold Dark Matter model to constrain the quasar 
halo masses. The same can also be done for high redshift 
galaxies where the quasar sightlines pass close to galaxies. 
Here the proximity effect is known definitively to be too 
small to affect the spectra (see e.g. Bruscoli et al. 2003) and 
instead increased absorption is also seen (Adelberger et al. 
2003). Some related work on the Lya optical depth around 
quasars has been carried out by Rollinde et al. 2005 and the 
distribution of damped Lya absorption around foreground 
quasars by Hennawi et al. (2006). 

Looking at the profile of Lya forest absorption around 
quasars in two dimensions, one parallel and one perpen- 
dicular to the line of sight allows one to quantify redshift 
distortions. As the gravitating mass governs the amount of 
squashing seen (e.g., Kaiser 1987, Regos & Geller 1989), 
measuring it may allow us to infer the mass enclosed within 
a given radius in a complementary way to looking at the 
mean strength of Lya absorption. 

Our plan for the paper is as follows: In §2, we describe 
the hydrodynamic simulations we used, as well as how Lya 
spectra were made. In this section we also compute the mean 
density and velocity profiles as a function of halo mass, di- 
rectly from the simulation mass and velocity data. We also 
compute the mean Lya forest flux as a function of quasar- 
pixel distance, for different mass bins. In §3 we describe 
how we fit the Lya forest flux-distance trend with a simple 
power-law model and how the power law parameters depend 
on the halo mass. In §4 we describe the data from the SDSS 
and the LBG data from Adelberger and use our results from 
§3 to constrain the halo masses. In §5, we examine the red- 
shift space anisotropy of clustering. In §6 we discuss our 
results and conclude. 



2 SIMULATED QSO SPECTRA 
2.1 Simulations 

We use two large N-body + hydrodynamics simulations of a 
A CDM cosmology to make our spectra. The two simulations 
were run using the code GADGET-2 (Springel, Yoshida & 
White 2001, Springel 2005), and are described more fuly in 



Nagamine et al. 2005. The simulations have the same cos- 
mological parameters ( Qa — 0.7, fi m = 0.3, fib = 0.04 and 
Ho = 70 km s -1 Mpc -1 , <T8=0.9) but have different box 
sizes and particle numbers. This means that they have sig- 
nificantly different mass and spatial resolutions so that we 
can use them for a resolution study, in order to make sure 
that our results have converged. Both simulations include 
gas physics, heating, cooling, a prescription for star forma- 
tion (Springel & Hernquist 2003) and the effect of stellar 
feedback on has properties. 

One simulation is the G6 run, a cubic box of 100 
ft _1 Mpc on each side, with 486 3 dark matter particles and 
486 3 gas particles. This results in an initial mass per gas par- 
ticle of 9.7 x 10 7 M Q and dark matter mass 6.3 x 1O S M . The 
other is the D5 run, a cubic box of side length 33.75 /i _1 Mpc 
and 2 x 324 3 particles. 

2.2 Mock spectra 

We use simulation outputs from redshift 3 to make mock 
Lya spectra. This was done in the usual manner, by inte- 
grating through the SPH kernels of the particles to obtain 
the neutral hydrogen density field, and then convolving with 
the line of sight velocity field (see e.g., Hernquist et al. 1996) 
The mean UV background radiation intensity has been nor- 
malized so that the spectra have a mean effective optical 
length r e ff = 0.4 (where r cff = — ln(F)) or equivalently a 
mean flux (F)=0.67. We make 5000 spectra from each sim- 
ulation, where the sightlines are parallel to the box axes and 
the x — y positions of the sightlines are picked randomly. 

2.3 Quasar hosts 

For halo finding, we use the friends-of-friends (fof) method 
(e.g., Huchra & Geller 1982). A particle is defined to belong 
to a group if it is within some linking length (b) of any 
other particle in the group. We select clusters by using b = 
0.2 where b is the linking length as a fraction of the mean 
particle separation. This definition of dark matter haloes 
was shown by Jenkins et al. (2001) to yield a mass function 
with a univeral form. We use these dark matter halos as our 
quasar hosts in the simulation. We bin the quasar hosts in 
terms of their mass. For each simulation we use has 5 mass 
bins, logarithmically spaced. There are 4 overlapping bins. 
D5 (G6) has one more bin smaller (larger) than the other 
4. This is because in the D5 simulation there are too few 
halos with masses > 1O 12 M0 and in the G6 simulation the 
mass per particle is too large to resolve halos with masses 
< 1O 9 M . 

2.4 Mean baryonic density and infall velocity as a 
function of quasar-pixel distance 

The baryonic density around quasars is related to the 
amount of neutral hydrogen which absorbs light via the Lya 
energy transition. As we have mentioned in §1, baryons fall 
into gravitational wells formed by dark matter, and so the 
baryonic distribution will mimic that of dark matter. With 
two quasars with close angular positions, but at different 
redshifts, we can therefore investigate the distribution of 
dark matter by studying the Lya absorption lines around a 
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Figure 1. Baryonic density in units of the cosmic mean as a 
function of r (distance from a dark matter halo) in real space for 
the D5 simulation(top) and G6 simulation (bottom). We show 
results for different dark matter halo masses. 



foreground quasar. Quasars have high luminosities despite 
their high redshift, and so allow us to probe the high redshift 
Universe. 

Because the density profile around quasars in the sim- 
ulation is directly available to us from the simulation data, 
we examine this first, before moving on to study the Lya 
absorption profile. In order to compute the baryonic density 
profile, we loop over all quasars in a mass bin and all pixels 
in the spectra. For each quasar-pixel pair we compute the 
separation, r in comoving ft -1 Mpc. We average the density 
values for all pairs at a given separation and show the results 
in Figure [T] 

We naturally expect that the baryonic density should 
be a decreasing function of r and that the overall amplitude 
should monotonically increase as the quasar host mass in- 
creases (Figure [1]) We find that this is generally the case, 
except for the fact that the the smallest mass bin in each 
simulation shows a higher density profile amplitude than ex- 
pected (we will see later that this also the case for the mean 
flux profile) . When we look at the overlap region for the two 
simulations we will see later that this is clearly a resolu- 
tion effect. For the bins in quasar host mass which are well 
resolved in the two simulations there is good agreement. 

We also compute the mean infall velocity as a function 
of quasar-pixel separation. The infall velocity (i>i n faii)is de- 
fined as the velocity component pointing towards the center 
of a quasar host. It should depend on the gravitational pull 
of the mass surrounding a quasar and therefore decrease as 
r does and increase as quasar host halo mass increases (Fig- 
ure 0. The tangential component will add up to zero if we 
average over the velocity of gas around a quasar. As seen in 
the case of the baryonic density, the resolution effect may 
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Figure 2. Infall velocity as a function of r in real space for the 
simulations D5 (left) and G6 (right). Note that the plot scales 
are different because of different simulation size. The i>i n f a ii curve 
for the largest mass bin (M) = 2.9 X 10 10 in the D5 simulation is 
shown as bold lines with crosses. 

be responsible for some mass bins having higher Vinfall than 
predicted. There is discrepancy between the two runs for the 
same mass bin on large r scales. This is likely to be because 
the box size of the D5 run is much smaller and peculiar 
velocities are sensitive to large-scale density fluctuations. 

2.5 Mean Lya forest flux as a function of 
quasar-pixel distance 

For each mass bin, we calculate the mean absorbed flux de- 
fined as: 

<F(r))=]Tcxp-^\ (1) 

where r(r) is the optical depth for Lya absorption in red- 
shift space, in a pixel at comoving distance from the quasar 
r. The observed flux is the fraction of photons with a given 
wavelength that are left unabsorbed by neutral hydrogen. 
Since the Lya absorption should reflect the presence of neu- 
tral hydrogen, we expect the observed flux to be correlated 
with the baryonic density, which increases close to quasars 
(small quasar-pixel separation r) as shown in the previous 
section. 

We use will use the G6 run, averaging over 500 spec- 
tra in this section as its mass resolution is suitable for the 
halo mass range relevant for quasars. We plot the mean flux 
in Figure [3] where we see that the flux asymptotically ap- 
proaches a mean value on large scales, which is 0.67 in this 
case, the value set in the simulation. The bottom panel of 
Figure |3] shows the case when we add a Gaussian redshift er- 
ror of 150 km s _1 to the quasar position. This is to reflect 
the measurement errors on observational determination of 
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Figure 3. The mean transmitted flux in the Lyct forest in red- 
shift space, as a function of quasar-pixel distance. Top: no red- 
shift error has been added. Bottom: a Gaussian redshift error of 
150 km s _1 has been added. The flux asymptotically approaches 
(F)=0.67, the mean value that has been set in the simulation. 



the quasar redshift. We can see that in this case, the shape 
of F(r) is different on small scales, r < 1 h~ 1 Mpc and also 
scatters more, especially for smaller mass halos. Since this 
is the G6 run, it may not be accurately reflecting the be- 
haviour of quasar hosts of small mass in any case. On larger 
scales, which will be relevant for our fitting later, there is 
not much change. 

The plot shows that the overall flux level F(r) decreases 
with mass, because more baryonic matter near QSOs means 
more absorptions of photons by the Lya forest, which re- 
sults in less observed flux. As with the baryonic density (Fig- 
ure [T]) , we see that the smallest mean flux is not from the 
smallest mass bin. However, we expect this to be due to res- 
olution effects. Apart from this, there is a good correlation 
between absorption and halo mass. 



2.6 Mean flux as a function of quasar-pixel 

separation perpendicular and parallel to the 
line of sight 

The flux in redshift space is subject to redshift distortions 
because of peculiar velocities, which act in a way similar to 
their effect on the autocorrelation function (Kaiser f987). 
The baryons around a quasar will infall towards it because 
of the gravity from the overdense region. This will result 
in the compression of the flux profile in the line of sight 
direction. The squashing effect will be greater for quasars 
with of larger host mass as they are associated with larger 
density fluctuations. We decompose the quasar-pixel separa- 
tion r into the transverse(cr) and line of sight (ir) directions 



it ) and plot the flux as a function of these two 
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Figure 4. The mean observed flux in terms of transeverse and 
line of sight directions. The flux goes through redshift distortions 
due to peculiar velocities and looks squashed in the LOS direction. 
From the top left: 5.1 xlO 9 M , 2.9x 10 10 M Q , 1.6xlO n M ,9.1x 
10 11 M(7> and 5.1 X 10 12 M P 
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mass is not very obvious, the distortion effect is clearly seen. 
In §5, we measure the anisotropy of the flux from simula- 
tions in redshift space by using Hamilton's formulae (1992) 
to quantify it and compare it with observational data. 



3 FIT TO THE OBSERVED FLUX VS 

QUASAR-PIXEL SEPARATION TREND 

In this section we parameterise F(r) using ro and 7 which 
are commonly used in the power-law auto correlation func- 
tion regularly applied to galaxy clustering data. Our aim is 
to eventually estimate the mass of SDSS quasar host ha- 
los and of Lyman Break Galaxy halos by using the relation 
between the fitted values of ro and halo masses. 

We use 5000 spectra from the G6 run, computing the 
mean flux as a function of r for 5 different mass bins. To 
fit the F(r) we need to use a covariance matrix, as the bins 
are correlated. In order to make a covariance matrix, we 
divide the simulation box into 125 subvolumes. We calculate 
a full covariance matrix (cr* 1 ) from the mass bin that has 
most quasars. For the rest of mass bins there are not enough 
quasars to compute a covariance matrix that is noiseless 
enough to invert. Instead we infer the full covariance matrix 
by scaling up from the matrix computed from large number 
of quasars, making use of the diagonal elements, which can 
be calculated. We assume that the relative cross-correlations 
are the same for each covariance matrix. We therefore only 
calculate the diagonal elements (er™) directly and calculate 
the off-diagonal elements from afj in the following way: 




parameters in Figure [4] Although its dependence on quasar 



(2) 



We use the mean flux as a function of quasar-pixel difference 
computed for 5 different mass bins. We then make curves for 
a total of 80 mass bins (18 bins in-between each) by linear 
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Figure 5. The mean flux as a function of quasar-pixel distance 
from the G6 simulation (symbols) and fits using equation|3](lines.) 



interpolation. We also build a covariance matrix accordingly 
for each sub-mass bin, so that we can compute a smooth 
trend of fit parameters with mass. 



3.1 F(r) fit as a function of halo mass 

In order to find an equation that best describes F(r), we try 
a power law form that is commonly used to fit the two point 
correlation function: 
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-Tea 1 + 



(3) 



This formula describes the mean flux as an increasing func- 
tion with r that asymptotically approaches a certain value 
(because there will be no relative enhancement of neutral hy- 
drogen far from a quasar). As r becomes large, the second 
term in the exponent will decrease, leaving the the flux close 
to exp _Teff which is the asymptotic value, and when r r is 
small, Lya absorption depends strongly on the distribution 
of matter around a quasar. 

The value of T e g is set so that the asymptotic value of 
equation [3] matches the one from the simulation. We fit the 
observed flux from the G6 simulation to equation|3]and show 
the results in Figure [S] These results show that even though 
it is very simple, equation [3] gives a very good description 
of the observed mean flux as a function of scale, at least for 
the r range of interest. 

We have fitted ro and 7 for the 80 mass bins, carrying 
out a \ 2 analysis using the full covariance matrix built as 
explained previously. Figure [S] shows the best fit values of ro 
and 7 for the D5 and G6 simulations as a function quasar 
mass. Both ro and 7 decrease at first and then go up, which 
is consistent with the non-monotonic behaviour of the mean 
flux in terms of mass as shown in Figure [3] As we mentioned 
before, this is due to resolution effects, as we can see by for 
example examining what happens to halos which approach 
the resolution limit of the G6 simulation, and compare them 
to the well resolved D5 halos of the same mass. 

We see that the slope is getting steeper when quasar 
mass increases, which implies that the overall amplitude of 
the flux on small scales decreases faster. Since the G6 run is 
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Figure 6. The best fit ro and 7 fitted to the quasar host-Lya 
forest profile vs. quasar host halo mass. 



more suited to the mass range of quasar hosts, we use the 
ro-mass trends from this run in our fitting , although our 
conclusions are not sensitive to this. 

We also carry out the same fitting procedure with the 
G6 run, after a Gaussian redshift error with a — 150 km s _1 
is added to the quasar redshifts. Both ro and 7 decrease with 
the addition of the redshift errors (also shown in Figure [6]). 
The ro value is not much affected by the redshift errors, 
however. 



3.2 Application to observational data 

Our next step is use the ro-mass relation to estimate the 
mass of dark matter halos for quasars observed in the SDSS, 
and LBGs, both at redshift z — 3. We will use equation [3] 
to find the ro and 7 values that best describe the mean flux 
as a function of separation and then use the best fit ro to 
estimate the mass. 



4 ESTIMATED MASS LIMIT FOR SDSS 
QUASARS 

We use ~3000 observed quasars from SDSS data release 3 
(Abazajian et al. 2005), the subset of the total data that 
have redshifts z > 2.4. This results in a mean redshift for the 
Lya forest pixels used in our analysis of z = 3.0, matched 
to the simulations we have been using. This quasar sample 
is rather bright, with mean G magnitude -27.5. It would be 
ideal to have many pairs of close-by quasars but unfortu- 
nately they are rare in observations. Of the 3000 quasars, 
only ~ 100 are close enough to contribute to the F(r) for 
r < 20 ft Mpc. In order to compute F(r), the pixels in 
the SDSS spectra we convert the pixel redshifts and quasar 
positions into comoving cartesian coordinates assuming a 
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A CDM cosmology and then compute quasar-pixel separa- 
tions. Only foreground quasars are used in our analysis, so 
that we are not sensitive to the regular proximity effect. The 
procedure is the same as in Croft (2004) , where it is outlined 
in more detail. 

The results for F(r) for the SDSS quasars are shown 
in Figure [7] along with the best fit power law curve and 
jackknife error bars. In order to carry out the fit, we use the 
covariance matrix derived from the simulations, scaled so 
that the diagonal elements are the same as the observational 
error bars. 

Using the ro vs. mass trend for the G6 simulation, we 
find the estimated mass for SDSS quasars host dark matter 
halos is log M/M Q =f 2. 34to.73. Although not used in the 
analysis directly, fitting 7 helps narrow down the quasar 
mass. We also tried fixing a 7 (not necessarily the best fit 
value) and varied only ro but could not recover a mass 
limit. In the case of the Gaussian error of 150 kms -1 , 
added to the quasar redshifts, the estimated mass limit be- 
comes logM/M Q = 12.48lE!;73. The result agrees well with 
that when the redshift error is not included, making the best 
fit mass slightly larger, by ~ 20%. Given the large error bars 
on the halo mass, this is not significant. 

The details of quasar formation and evolution in time 
are still in debate, which makes it difficult to predict the 
mass of quasar hosts, and there are few estimates their these 
masses from observational data. Wyithe & Padmanabhan 
(2006) report that mass estimates of quasar host dark mat- 
ter haloes are from I0 11 Mq — 1O 13 M0, depending on the 
time-dependence of their evolution. Our value lies within 
this range of results. 



4.1 Estimated Mass Limit for Lyman Break 
Galaxies 

We apply our method to Lyman Break Galaxy data (Adel- 
berger et al. 2003) . The reported mass of LBGs from other 
methods at z = 2.9 is 10 112 - 10 11 ' 8 (Adelberger et al. 
2005) while other groups have slight different estimates: 
10 ii.6±o.3 M@ ( Weatherley & Warren 2005) and 10 n ' 3 M Q 
(Somerville et al. 2001). 

The fitting procedure is same as previously used except 
that we ignore the the first three bins because the data for 
small r is not consistent with the trend in the simulation 
despite the large error bars(Figure|8}. If the data points are 
valid, it would need further investigation, and indeed these 
results may be indicative of starburst winds disturbing the 
Lycv forest on small scales. But for simplicity, and because 
winds from galaxies are unlikely to propagate far, we only 
use regions greater than 1.5 h~ Mpc in our analysis. We fit 
the LBG data to the equation [3] and then estimate the mass 
using the ro-mass relation found in §3.1. The estimated mass 
of LBGs we find is logM/M Q =11.13±^. The best fit value 
of mass is consistent with Adelberger's estimate. 



4.2 Dependence on cosmological parameters 

In carrying out our analysis and computing constraints on 
the halo masses of quasars and LBGs we have assumed that 
the ro- mass relation from the simulations applies to the 
observations. This is probably a valid assumption, as the 
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Figure 7. The mean transmitted flux as a function of quasar- 
pixel distance, for SDSS DR3 quasars (crosses with error bars) 
fitted with Eqn. \3\ (dotted lines). 
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Figure 8. The mean flux of Adelberger's data (crosses with error 
bars) fitted to Eqn. [3] (dotted lines) 



cosmological model used in our simulation is close to that 
found from for example the WMAP satellite data (Spergel 
et al. 2003, 2006). As the error bars on our best fit quasar 
halo mass are large, we expect errors on cosmological pa- 
rameters to have relatively little effect. In order to verify 
this, we vary the least constrained parameters, the ampli- 
tude of mass fluctuations, erg, running new simulations with 
different amplitudes and seeing if we get the same quasar 
host masses if we use these to compute the ro-mass relation. 

The new simulations are of dark matter only, and it 
is assumed that the gas traces the dark matter distribution. 
Three runs are made, with <T8=0.7, 0.9 and 1.1. They are run 
in a 50 ft -1 Mpc box, using 256 3 particles. Spectra are made 
from the dark matter distributions in the manner described 
by e.g. Croft et al. (1998). Although the simulations do not 
include any hydrodynamics we expect the relative values of 
ro for the different eg runs to closely approximate the values 
that would be obtained with full hydro simulations. 

The fitting procedure is the same as the D5/G6 sim- 
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Figure 9. ro and 7 vs. quasar mass for different ag values. (0.7, 
0.9, and 1.1) No hydrodynamics is included in the simulation. No 
redshift error has been added. 

ulation cases except that the covariance matrix was built 
from 27 subvolumes because the simulation size is smaller 
than the G6 run. We fit the quasar data from SDSS using 
equation [3] and find the mass limits using the newly found 
ro vs. mass relation. We use no redshift error in the analysis. 
We find the mass limits as follows: logM/M Q =12.21±°'g4, 
logM/M Q =12.88±°;^, and logM/M =12.68±°;^ for a 8 = 
0.7, 0.9, and 1.1 respectively. The comparison with the re- 
sults in §3.1 (logM/M Q = 12.34lo;73) indicates that the mass 
fluctuations do not have a significant effect on estimating 
quasar masses, at least within our error bars (less than a la 
effect). 

4.3 Dependence on r e « and temperature 

The hydrodynamic simulations have a gas temperature at 
the mean density close to 20,000 K (see observational de- 
terminations by e.g. Schaye et al. 2000) and mean effective 
optical depth (r c ff = — In {F))= 0.4, which has been set to be 
consistent with observational data. Both of these quantities 
have uncertainties, and these may affect our recovery of the 
quasar dark matter halo mass from Lya forest clustering. 
We use the dark matter only simulation to estimate what 
the total errors on halo mass will be by quantifying the effect 
of changes in these parameters. We make Lya spectra from 
the simulation setting To = 10, 000 K and 30,000 K instead 
of 20,000K. We also make spectra by changing the the mean 
flux (F) by ±10% from our original value of (F) = 0.67 In 
these cases the mean transmitted flux becomes 0.603 and 
0.737, which translates to mean effective optical depths T e s 
of 0.305 and 0.506, respectively. This range of T e g covers the 
values measured by T e g by Rauch et al. (1997) and other ob- 
served values: For example McDonald et al. (2000) reported 




Mass (M s h l ) 



Figure 10. ro and 7 vs. quasar mass when mean observed fluxes 
are 90 and 110 % of the nominal value (0.67) and the gas tem- 
perature is 10,000K and 30,000K. The solid bold lines represent 
(F(r)) =0.67 and T=20,000K case (<r 8 =0.9) No hydrodynamics 
is included in the simulation. No redshift error has been added. 



r off 0.380 ± 9% at z=3 and Schaye et al. (2003) obtained 
(F) = 0.696 at z=3. Although we vary (F) by ±10%, this 
range is much larger than the scatter in the observed val- 
ues, which we conservatively estimate to be ±5% at the 1 a 
level. We use this value when we compute the effect of (F) 
uncertainty on the quasar host halo mass. 

The fitting procedure is the same as in the previous 
sections. Figure [TOl shows the mass vs. ro relation for each 
of these cases. Adding the uncertainty in quadrature for 
asymmetrical error bars, we find: \/0.285 2 ± 0.41 2 ± 0.19 2 = 
0.53 for the upper error bar and V0.49 2 ±0.73 2 +0.11 2 = 
0.89 for the lower error bar, where each term represents 
errors from the mean flux, statistical errors (estimated 
using the G6 run), and temperature uncertainty respec- 
tively. Including the 150 kms' 1 velocity uncertainty, our 
mass estimate for z — 3 SDSS quasar host halos be- 
comes log M/M0=12.48^Q gg . This value is in good agree- 
ment with the reported mass of quasar host dark halo 
(Mdmh) Mdmh = 3.0 ± 1.6 x 10 12 /i _1 M© inferred by 
Croom et al. 2005 from the autocorrelation function of 
quasars in 2dF data. It is also within 1 a of Mdmh = 
5.2 ± 0.6 x 10 12 /i _1 Mpc recently obtained by Myers et al. 
(2006) from clustering of 300,000 photometrically classified 
SDSS QSO. In the same manner, when we add the (F) er- 
rors and temperature errors to the LBG mass estimate, we 
obtain log M/M Q = 1 1.13^0 55. 
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5 MEASURING THE REDSHIFT SPACE 
ANISOTROPY OF CLUSTERING 

5.1 Quadrupole to Monopole Ratio 

The galaxy autocorrelation is distorted in redshift space be- 
cause the line of sight direction is displaced due to pecu- 
liar velocities. The overall flow towards a galaxy mass en- 
hancement will result in a squashed correlation function. 
The anisotropy on large scales can be parametrized in lin- 
ear theory using (3 = /b where b is the bias parameter 
(Kaiser 1987, Hamilton 1992.) Hamilton (1992) derived the 
expected value of the quadrupole to monopole ratio of the 
correlation function in terms of /3, using linear theory. Even 
on non-linear scales, this ratio can be used as a measure 
of the squashing effect. This measure has been employed 
by both the SDSS and 2dF survey groups to quantify the 
anisotropy of galaxy clustering in redshift space, (e.g., Ze- 
havi et al. 2002, Hawkins et al. 2003) Replacing the correla- 
tion function £(r) in Hamilton (1992) with our quasar-Lya 
forest cross-correlation function, F(r) = F(r)/ {F) — l where 
(F) — 0.67, we will compute the quadrupole to monopole ra- 
tio and see if it has a dependence on quasar host halo mass. 

Following Hamilton's formulation, we decompose the 
mean flux F(r) into F(r,fi) using Legendre Polynomials: 



F(r,n)=J2 F i(r)Pl(»), 



(4) 



where p is the cosine of the angle between the line of sight 
and the pair separation vector in redshift space, and Pi are 
Legendre Polynomials. Integrating over all cosine angles, we 
get: 



Fi 



F(r,/i)(l + 2l))fl(/i)dAt, 



(•») 



We calculate the quadrupole moment of the flux as follows 
(Hamilton 1992): 



Q(s) = 



(4/3)0 + (4/7)/3 2 
l + (2/3)/3 + (l/5)/3 2 

Ms) 

(3/ S 3) f° Fo (s>)s>2ds>-F (s) 



(6) 



We carry this out for the G6 simulation and also for the 
SDSS quasar data. Our expectation is that the magnitude of 
quadrupole moment is greater for more massive quasar host 
halos. If this is the case, the goal is then to use the relation 
between Mdmh and the quadrupole moment to constrain 
the mass of SDSS quasars in a manner independent of the 
overall strength of absorption which we have use in the main 
part of this paper. 

In Figure[TTJwe show the quadrupole to monopole ratio 
for halos in several mass bins. Contrary to our expectation, 
the quadrupole moment does not show a very strong depen- 
dence on Mdmh (Figure ITT]) . It is not possible to see any 
useful monotonic relation between Q(s) and Mdmh- 

We have computed Q(s) from our SDSS quasar sample, 
and find results that are consistent with the simulations, but 
with extremely large error bars. From this study, it seems as 
though the simulations reproduce the anisotropy of cluster- 
ing, but there is not any expectation of being able to use its 
magnitude to measure quasar host halo mass. On the other 
hand, the anisotropy will be sensitive to assumed cosmology, 




Figure 11. Quadrupole moments of the mean flux from the G6 
simulation and SDSS data (stars with error bars). No redshift 
error has been added. 



and a much larger sample of this type of data could be used 
to carry out an Alcock-Paczyhski (1979) test. 



6 SUMMARY AND DISCUSSION 
6.1 Summary 

In this paper, we have explored the profile of Lya absorption 
with distance around quasars using pairs of close-by quasars 
in high resolution hydrodynamic simulations. We have found 
the following: 

• The mean Lya forest transmitted flux from a back- 
ground quasar passing through a region close to a foreground 
quasar and the infall velocity of matter around a quasar are 
both dependent on quasar mass. 

• We can fit the observed flux vs quasar-pixel separation 
r to an equation exp(— r a g (l + (^) -7 )) by varying ro and 7. 
This simple power law form, derived from fits to the galaxy 
autocorrelation function works very well. 

• There is a strong dependence of ro in this fit to the 
mass of the quasar host halo in the simulation. From this ro 
vs. mass relation, we can estimate quasar halo mass. 

• The estimated mass of SDSS quasars at z — 3 
is logM/M =12.48t° gg. We do not see any significant 
changes in the result when we include redshift errors of 
150 km s _1 or a different amplitude of mass fluctuations. 

• The estimated halo mass of Lyman Break Galaxies can 
be found using the same technique. Fitting to the data of 



Adelberger (2003) we find logM/M = 11.13±° ™. 

• The squashing of flux in redshift space can be quan- 
tified using the quadrupole to monopole ratio. There is a 
no measureable dependence of the squashing effect on halo 
mass, although the anisotropy seen is consistent with obser- 
vational data. 



6.2 Discussion 

Due to their great distance from us, it is difficult to mea- 
sure the dark matter halo masses of high redshift quasars. 
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For example, few background galaxies exist to make lens- 
ing measurements possible, and the rotation curves of their 
host galaxies may require larger telescopes than currently 
available to capture. 

The strong level of Lya forest absorption around bright 
z — 3 SDSS quasars that we have seen indictates that they 
have halo masses ~ 3 x 1O 12 M0, comparable to elliptical 
galaxies at the present. This is not unreasonable because 
these bright quasars are much rarer than present day el- 
lipticals with this mass. Our mass estimate is roughly 20 
times larger than that for the dark matter mass of LBGs at 
the same redshift, indicating that it is unlikely that bright 
quasars randomly sample the population of LBGs as hosts. 

In our analysis, we have ignored the transverse prox- 
imity effect, the possiblity that foreground quasars will de- 
crease the amount of Lya absorption with local photoion- 
ization from their radiation. If this effect is present, it is 
certainly not dominant, as we have seen that instead there 
is more absorption close to foreground quasars. If it does ex- 
ist, but is subdominant, then the true amount of absorption 
without this extra radiation should be greater and we are 
underestimating it and hence the halo masses. ~ 3 x 1O 12 M 
therefore represents a lower limit to the dark matter mass 
of quasar hosts. There are however several ways, including 
quasar beaming and short lifetimes mentioned in §1 that 
the transverse proximity effect might not be relevant in any 
case. 

The redshift distortion of the absorption is unfortu- 
nately not promising as a complimentary method for con- 
straining the dark matter halo mass. However the insensitity 
of the distortion to mass means that the correlation function 
it could be a useful probe of cosmic geometry. As the Lya 
forest- quasar clustering signal appears to be much stronger 
Lya forest two point clustering, this may be an efficient 
way of carrying out a high-z Alcock Paczyriski test. 
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